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Abstract 

A Monte Carlo study of two jet observables that can be used to obtain 
jet samples enhanced in quark or gluon content is presented. We compare 
the performance of one of these methods on pure Monte Carlo events 
with that after the transport of events through the full simulation of the 
ALICE experiment. It is shown that the performance of the method is 
■ not deteriorated by experimental effects like resolution and inefficiency. 

1 Introduction 

. The first step towards the use of jets as probes of the properties of the hot and 

' dense medium produced in heavy-ion collisions is to characterize the properties 

QO ■ of the jets produced in vacuum. This characterization is by itself an important 

lO ' test of perturbative QCD and will also help to tune Monte Carlo generators at 

. LHC energies. 

We are interested in studying the response of the medium to quark and 
gluon jets independently in order to test the theoretically predicted differences 
in energy loss [T] [5] due to their different color charge factor. The objective 
of this study is to test the feasibility of different methods, in proton-proton 
collisions, that allow to obtain jet samples enhanced in quark or gluon content. 



2 ALICE charged jet reconstruction 

ALICE[5][1][5] is a dedicated experiment at the LHC to study heavy- ion colli- 
sions. It has excellent tracking and particle identification capabilities j3], that 
makes it an ideal tool to study jet properties like fragmentation. 

The jet reconstruction presented here is based only on charged particles 
within the ALICE acceptance \r)\ < 0.9. In the future, ALICE will also be able 
to perform full jet reconstruction with the combination of tracking from the 
central barrel detectors and calorimetry with the EMCal detector. 



* Preprint of an article submitted for consideration in IJMPE ©2010 [copyright World 
Scientific Publishing Company] http://www.worldscinet.com/ijmpe/ijmpe.shtml 
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3 Data sample and jet finder 



The data sample used for this work consists of jet events in proton-proton colli- 
sions at -y/s = 7 TeV produced with Pythia[6 using the set of parton distribution 
functions CTEQ4L 7,. Full detector simulation of these events was performed 
using AliRoot[3] in order to estimate the performance of the methods under 
realistic conditions. 

The jet finder used was the UA1 cone algorithm, using a cone radius R = 
y Arj 2 + A0 2 of 0.4, a seed of 4 GeV and requiring a minimum jet transverse 
energy E^ C U of 5 GeV. In the following, the parton association to jets was done 
by selecting the most energetic parton inside a subcone with R =0.3 around the 
reconstructed jet axis. The quark/gluon content of the sample is determined 
this way. The charged jet reconstruction was found to introduce a bias towards 
finding a larger fraction of the quark initiated jets compared to the full jet 
reconstruction. This bias was found to increase due to the use of a small cone 
radius. 



4 Quark and gluon enhancement methods 

In the following subsections, the two enhancement methods studied will be 
described followed by a comparison of their performances. 

4.1 Second central moment of transverse structure 

The second central moment of the transverse structure [8] is defined as in Eq. (fTJ. 
For each component a (a = <j>, rj) the second central moment is obtained using 
Eq. ([2} and Eq. ([3]), where the index j corresponds to the rj and <p of the jet axis 
and the index i runs over the charged tracks that are part of the jet, having a 
transverse momentum px ■ 

(5Rl) = (64) + (8vl) (1) 
(Saf) = <^ 2 ) " (8a,) 2 (2) 

(8a-) = S ^~ Q f x Pt (3) 

Figure [T] shows the distribution of the mean values of 8R 2 as a function of 
E^ ct for quark, gluon and all jets. The tracks used to perform this analysis were 
charged tracks with a pt > 1 GeV/c and with \rj\ < 0.9. As one can observe 
from the distributions, gluon jets have on average a larger mean value of 8R 2 
than quark jets over the whole range of E^ ct studied. This difference in (SR 2 C ) 
is the basis of this enhancement method. 

lr The UA1 jet finder algorithm produces massless jets so p^, ct =E^, ct . 
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Figure 1: Mean values of SR^ as a 
function of E^ ct for gluon jets in red 
triangles, quark jets in blue squares 
and all jets in black circles. 



Figure 2: Mean values of NT90 as 
a function of E^ ct using the same 
symbols and colours as in Fig. 1. 



4.2 Track counting method 

The second method studied is based on a modified version of that proposed 
to be used with segmented calorimeters [5]. The track counting algorithm is 
implemented in the following way: the tracks that compose a given jet, in this 
case tracks inside the cone radius around the jet axis, are ordered in decreasing 
magnitude of transverse momentum. Then the transverse momentum of the 
ordered tracks is added until a certain fraction of E^ et is recovered, in this case 
90%. This minimum number of tracks is defined as the variable NT90. 

Figure [5] shows the distributions of the mean values of NT90 for quark, gluon 
and all jets. One can observe that the mean values of this variable are larger 
for gluon jets compared to quark jets of the same E^ ct . This way it is possible 
to select quark or gluon jets based on the multiplicity needed to recover 90% of 
the E^ ct . 

4.3 Comparison of results 

For both methods, the strategy to select quark jets was based on setting an 
upper value on the variables so the resulting samples would be enriched in 
quark content, and gluon jets were selected by requiring a minimum value of 
AT90 or (8R 2 C ). 

In order to compare the performance of both methods, their results for 
selecting quark jets at a given jet energy and achieving the same purity of 
the order of 60% to select quark jets were compared. The method of the second 
moment of the transverse structure is 30% less efficient than the track counting 
method, as shown in Table 1. 



Table 1: Comparison of discriminating methods for E^ et = 20 ± 5 GeV. 



Method 


Purity 


Efficiency 




61% 


38% 


NT90 


64% 


69% 
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In the next step of the study, we tried to combine cuts using both variables to 
increase the discriminating power. The result, however, is that both variables 
were highly correlated. In the following, only the performance of the track 
counting method will be described comparing results using only Monte Carlo 
(MC) events and the results after the full reconstruction chain. 

5 7VT90 method performance 

In order to evaluate the enhancement power of the variable NT90, the efficiency 
(e) and purity (P) for selecting quark/gluon jets was studied as a function of 
the cuts on NT90. This study was done for a fixed value of the E^ ot , the value 
was chosen to be 20 ± 5 GeV. 

Figures |3] and 0] show the results from the pure MC events and for the same 
events after the full detector simulation. 
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Figure 3: e and P for gluon jets as 
a function of the cut on AT90. The 
results from pure MC are shown 
with black lines and open triangles 
and the results after the ALICE 
simulation with a red line and full 
triangles. 
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Figure 4: e and P for quark jets as 
a function of the cut on NT90. Re- 
sults from pure MC are shown with 
black lines and open triangles and 
the results after the ALICE simula- 
tion with a blue line and full trian- 
gles. 



For gluons jets in Figure |3l the first bin corresponds to not performing any 
cut. By setting the cuts, the minimum value on the number of tracks necessary 
to recover 90% of the E^ ct is fixed. For the quark case, shown in Figure |4j the 
last bin represents the result when no cut is applied and successive cuts are 
made towards the first bin, in this case setting an upper value for the value of 
NT90. 

As one can see, the effect of the detector response does not affect the per- 
formance of the method. Table 2 shows an example of the purity achieved by 
setting cuts on NT90, but one should keep in mind that the method is better 
suited to disentangle jets produced by quark or gluon with higher E^ ct . It is 
in the high energy regime where the differences in the properties of quark and 
gluon jets are larger using NT90 (see Figure [2]) . 
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Table 2: Purity and efficiency for quark and gluon samples for E^ ct = 20 ± 5 



GeV. 



Parton 


Cut value 


Purity (enhancement) 


Efficiency 


Gluon 


7 


90% (14%) 


18% 


Quark 


4 


36% (12%) 


63% 



6 Conclusions 

The feasibility of a method to obtain samples enriched on quark or gluon content 
has been studied. It was shown using a full simulation of ALICE that it is 
possible to use a tagging variable based on charged track multiplicity to obtain 
jet samples enhanced in quark or gluon content in proton-proton collisions. For 
20 GeV jets it is possible to increase the purity of gluon jets up to 90% and of 
quarks jets up to 36% having efficiencies of 18% and 63% respectively. Using 
higher E^ ct it would be possible to obtain samples with higher purities of quark 
jet content. 

It is clear that the 7VT90 tagging method biases the jet samples towards 
having a certain structure, i.e. soft gluon jets and hard quark jets. However 
iVT90 can also be used as a jet structure observable that can be studied to 
compare jets from proton-proton collisions with those produced from heavy-ion 
collisions. 
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